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Abstract 

We study a class of three-point functions on the de Sitter universe and on the asymptotic 
cone. A blending of geometrical ideas and analytic methods is used to compute some remarkable 
integrals, on the basis of a generalized star-triangle identity living on the cone and on the complex 
de Sitter manifold. We discuss an application of the general results to the study of the stability 
of scalar particles on the Sitter universe. 



1 Prologue 

The main result of this paper is the following formula: 

h d ( K , „, a) ^ J~ pzrfju) p:rfju)p:ffju)(u* - 1)-^ d U = (i) 

d TT "P / d—1 i ieK-\-ie' v-\-ie" \ \ 
22 iL,e',e"=±l l {— H 2 J 

~ (47r)§r [n e=±1 f (E± + ien)] [U^ ±1 f + ^ f ( O + ie»\)] (2) 

expressing the integral of a product of three Legendre functions of the first kind [T] as a ratio of 
products of Euler Gamma functions. This beautifully symmetric formula is not listed in any of 
the handbooks on integrals of special functions available to us and appears to be new. The steps 
involved in the proof-computation also give rise to many interesting quantities having possibly ge- 
ometrical interpretations that we have not yet fully explored in their mathematical and physical 
consequences. 

The problem of computing the integral ([I]) originates from the study [2j [3] of particle decays in a 
de Sitter universe. In that context, the real parameters n, v and A are related to the masses of the 
particles involved in the decay process and d is the (complex) dimension of the de Sitter universe 
where the process takes place. However, the methods and the results we are going to present have 
presumably a wider interest and range of applications. 
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The study of particle decays in a de Sitter universe was initiated by 0. Nachtmann [4] in 1968. He 
showed, in a very special case, that while a Minkowskian particle can never decay into heavier products, 
a de Sitter particle can, although this effect is exponentially small in the de Sitter radius. 

The subject has acquired a greater physical interest with the advent of inflationary cosmology. In 
particular, the idea that particle decays during the (quasi-)de Sitter phase may have important con- 
sequences on the physics of the early universe has been suggested recently [6] (see also the related 
works [3 [8]). The mathematical and physical difficulties related to the lack of time-translation sym- 
metry of the de Sitter universe, and more generally of non-static (cosmological) backgrounds, have 
been tackled in [51 E] by using the Schwinger-Keldysh formalism, which is suitable for studying cer- 
tain aspects of quantum physics of systems out of equilibrium. However, the approach described in 
[5j [6] necessitates the introduction of a practical notion of lifetime of an unstable particle which is 
completely different from the definition commonly used in quantum physics. 

Actually the lack of a commutative symmetry group of spacetime translations renders the mathematics 
extremely complicated but docs not prevent computing the inverse lifetime of an unstable de Sitter 
particle according to the usual definition, namely as the inclusive transition probability per unit time 
from an initial state to every possible final state: this computation, initiated by Nachtmann long 
ago, has been performed at first order in perturbation theory in [21 Ej . After writing the relevant 
perturbative amplitude, the computation of the lifetime of the de Sitter unstable particle amounts to 
two essentially distinct and independent steps: 

1. taking the adiabatic limit of the infrared regularized inclusive amplitude, i.e. removing the 
infrared cutoff coupling factor necessary to make the integral expressing the amplitutde converge; 

2. computing the so called "phase space" coefficient, a quantity which only depends on the masses 
of the particles involved in the decay process. 

The first step has been largely discussed in [51 [3] and the resulting mathematical structures elucidated 
there. The second step reduces to computing the integral at thr RHS of (TT|). 

The quantity h^K, v, v), relative to a decay into two identical particles, has been computed in [5J [3]. 
This special case already exhibits some concrete mathematical difficulties, and has been solved in a 
purely analytical way by the use of Mellin transform techniques and the evaluation of a Barnes-type 
integral [1[3]. 

The above method fails however to provide a solution for the general case of the production of two 
non-identical particles, i.e. fails to give a solution to the general integral ([1]). For odd values of 
d the Legendre functions of the first kind reduce to trigonometric-type functions; in these cases a 
direct computation of hd{n, v, A) is possible. To give an example, one can solve the three-dimensional 
problem (d = 3) by an elementary computation. Indeed for d = 3 there holds the particularly simple 
expression of the Legendre function: 



PJ. + J^v) = J — . (3) 

2 +lu V 7r sinn v v 

A straightforward computation then gives: 

1 sinh(7Tft) sinh(7rA) sinh(7rz/) 

] = vWA ^h^^ cosh lia±^l C osh 3^±3 cosh iia±£±!d ■ (4) 

The general odd-dimensional case d — 2n + 1 can similarly be tackled by (increasingly cumbersome) 
elementary integration techniques. On the contrary, the computation of the even dimensional cases 
(including the physically relevant four-dimensional de Sitter universe) is very very far from obvi- 
ous. 
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Some of the geometric ideas necessary to overcome the difficulties of the integral (TTJ were contained 
in an unpublished work by one of uqj [5] . Combining those ideas with a the geometrical properties of 
the complex de Sitter manifold provides a way to solve the problem. The result is displayed in Eq. 
©• 

Beyond the study of the of de Sitter particle decays, there are other potential applications of the 
formula ([2]) and of the methods used to derive it which include the study of tensor product of repre- 
sentations of non-compact groups, many new integral relations involving products of hypergeometric 
functions, other applications to de Sitter and/or anti de Sitter QFT etc.. 

2 Legendre functions and de Sitter Klein-Gordon fields: a 
short review 

The computation of the integral ([T]) requires several steps in which the geometrical features of the 
complex de Sitter manifold enter in a crucial way. The first important step consists in returning to the 
meaning of the Legendre functions of the first kind P(u) as two-point functions of quantum fields on 
a complexified de Sitter spacetime [10l [HI [T2l [13] . The variable u, appearing as integration variable 
in the r.h.s. of ([2]), is understood as a geometrical invariant u = u(x,x') relating two points x,x' of 
a (complex) de Sitter hyperboloid. This idea allows in particular a natural way for understanding 
many of the mathematical properties of the Legendre functions and gives also a simple procedure to 
build many of their integral representations. Here follows a short account of the construction. 

Consider a (d + l)-dimensional Minkowski spacetime M d+ i; an event x is parameterized by a set 
of inertial coordinates x°, . . .x d ; the scalar product of two events of M d+ i is the Lorentz-invariant 
product x ■ x' = x°x'° — x 1 ^ 1 — ... — x d x' d . The rf-dimensional de Sitter spacetime is represented as 
the one-sheeted hyperboloid 

X d = {x e M d+ i :x-x = x 2 = -R 2 } (6) 

embedded in M d +i- The Lorentzian geometry of the de Sitter manifold is induced by the causal 
structure of the ambient spacetime: 

V+ = U € M d+1 : i 2 = i ■ i > 0, e > 0}, (7) 
C+ = U € M d+ i : e = £ • £ = 0, £° > 0}; (8) 

the future cone V + of the ambient spacetime induces the Lorentzian global causal ordering on the 
de Sitter universe: a; is in the future of x' if and only if x — x' G V + . The forward light-cone C + of 

lr The problem studied in [9] was the construction of a generating function for the Clebsch-Gordan coefficients of 
the group 5(7(2) following a method originally due to Wigner. Given three spinors ipi = i^i^Vi) * = 1, 2, 3, one can 
construct the unimodular invariants [4>i4>j\ = £,iVj ~~ £jVi\ * ne quantity I = [i/'i^]" 3 \^2^>3\ ai H'3' l Pl] a2 is an SU(2)- 
invariant polynomial trilinear in the vectors of the representation Dj i = §i* +m, 7i* m + m) ! (ji — m) ! = Xi(m), 
where a\ +ci2 = 2j3, etc.. The coefficients of the above polynomial are proportional to the Clebsch-Gordan coefficients. 
To determine the constant C it is sufficient to compute the following average on the unit sphere 

J = f |J"| 2 (ffii(ffi2df23 =C 2 ( Jl J2 J3 } =C 2 (5) 
J 1 1 \ m 1 m 2 m :i J w 

with the parametrization § = cos j e 1 ? and rj = sin ^ e -1 '?. Since llV'iV'j]! 2 = sm2 = ( ^a** ) one can recognize 
here the triangular invariant introduced in Section (|5j and computed in Section H79I I. Here the calculation of the 
invariant J is done for non-integer values of the exponents (non-compact case). The method we will use generalize the 
partial integrations that were enough to solve the integer case in [9] by using the fractional calculus. 
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the ambient spacetime also plays the role of the space of momentum directions in de Sitter momen- 
tum space [121 114) . The de Sitter invariance group is the Lorentz group of the ambient spacetime 
50(1, d). 

A de Sitter generalized free field <fi is fully characterized by its two-point vacuum expectation value 
W(x, x') which is assumed to be be a local and de Sitter invariant distribution. Since there is no 
global de Sitter energy operator, a true spectral condition does not exist in the de Sitter spacetime; 
there is however a suitable replacement that can be formulated p~2] by moving to the complex de 
Sitter manifold 

*f ={zeMjfl:z-z = -R 2 } (9) 
and requiring that W(x,x') be the boundary value on the reals of a de Sitter invariant function 
W m (z,z') holomorphic in the tubular domain T~ x T + with slow increase properties at infinity, 
where 

T± = {(M d+1 ± iV+) n X d c) } = {z = x + iye X d c \ yy>0, sign^ ) = ±}. (10) 

de Sitter invariance can then be used to show that W{z, z') is actually maximally analytic, i.e. it is 
analytic in the domain 

A = {(z,z')eA^ c) xA^ c) : (z _ z ') 2 ^R + }. (11) 

For a thermodynamical interpretation of the above analyticity property, see |12[ I13| 1 1 5 j . By intro- 
ducing the de Sitter invariant variable 

z • z' 

£ = — — ( = u when real and greater than one) (12) 
Rr 

there holds a simple description of W in terms of a function w of the single variable £, namely 
w(C) = W(z, z'), holomorphic in image of the domain A 

n = {cec, C^-l-R+}. (13) 

For a Klein-Gordon field cj> with mass m > the two-point function must also be a bisolution of the 
Klein-Gordon equation 

{n x + m 2 )W m {x,x') = (□*,+ m 2 )W m (x,x') = (14) 

where □ is the Laplace-Beltrami operator relative to the de Sitter geometry. It is useful to introduce 
a dimensionless parameter v related to the mass m as follows: 

m 2 R 2 = +v 2 . (15) 

By abuse of language we will call v a mass parameter even if it is dimensionless. Given a complex v 
the corresponding two-point function W m {z, z') = W v {z, z') = w„(C) is written in terms of Legendre 
functions of the first kind as follows: 

WJz,z') = — V-^- 2 ! -{C 2 -I)-— P ^ t . 16 

V ; 2(2^)ii? d - 2 VS ^ -h+™^ J V > 

z, z' are events belonging to A; the normalization ensures that the canonical commutation relations 
hold with the correct coefficient. 

The range m > m c = (d — l)/2R corresponds to the principal series of unitary irreducible represen- 
tations of the de Sitter group (y real) while < m < m c corresponds to the complementary series (y 
imaginary). These restrictions ensure that the boundary value W m is positive definite and therefore 
a quantum theoretical interpretation is available. Note also the symmetry property 

W»(x,x')=W- v (x,x') (implied by PlffjO = Pl^TJC)) , (17) 
which holds for all v and will play a role in one structural aspect of the derivation of formula (fl"|). 
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2.1 Plane waves expansion of Legendre functions 

There exists a Fourier-type representation of the two-point functions (|16p which is of fundamental im- 
portance to understand the above properties, and generally speaking to understand de Sitter QFT. It 
is constructed by using a natural basis of plane- wave solutions i/v of the Klein-Gordon equation 



a 



d- 1 



2R 



R 



{KG) V M*) = 

which are parameterized by the choice of a lightlike vector £ £ C + as follows: 

d—l , .... 



(18) 



(19) 



These waves are well-defined and analytic in each of the tubes T + and T~. Then, for z S T~ and 
z' G T + the following Fourier-type (i.e. momentum space) representation of the two-point function 
holds true: 



W u (z,z') 



r(^§i + w)r(^i - iv)e-™ 



«(0> 



(20) 



In standard coordinates the (d — l)-form is written 



3=1 



(21) 



In (f2"0"l) . the (d — 1)— form under the integration sign is closed and 7 denotes any (d — l)-cycle in the 
forward light-cone C + which is homologous to the following cycle 70 : the support of 70 is represented 
as follows by the unit sphere, S<j-i (equipped with its canonical orientation): 



supp7o = S d _i = C+ n U : e° = 1} = U G ^ + = f 1 ' 



r = 1}. 



(22) 



With this choice coincides with the rotation invariant measure on §d-i normalized as fol- 

lows: 

= / «(0 - (23) 



r(l) 



2.2 Lobatchevski space and a remarkable representation of Legendre func- 
tions 

A specially important parametrization of the Fourier-type representation is obtained by evaluating 
([2"0")) at the purely imaginary events z = — iy G T~ and z = + iy' G T+; t/ and y' can be 
visualized as points belonging to a Lobatchevski space, modeled as the upper sheet of a two-sheeted 
hyperboloid: 

Hd = feel w : y 2 =yy = R 2 , y°>0}. (24) 
We will make use of the following spherical parametrization of H^: 

y(u, n) = R(u, n 1 \Jv? -I,..., nV" 2 - 1) (25) 

where it > 1 and n G §d_i; in these coordinates the Lorentz-invariant measure dy is written 

dy = R d {u 2 -l) ± ^dudn (26) 
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where dn denotes the rotation- invariant measure on the sphere E>d-i normalized as in Eq. (|23[) . With 
the above specifications, Eq. (J27J) allow us to write: 

Wv{-iy,iy) = 2 d + i TT d I (y-0 (Z-y) d M 7 (0 = 



2(2tt)1 y ' 



d-2 



Y 4 p :^Jyy')- ( 2 7) 



Here and in the following we have set R = 1; by choosing in particular 7 = 70 and j/' = (1,0,..., 0) 
so that y ■ y' = y = u > 1, we then get the following integral representation: 



d-2 d-2 I 

/ 2 — ~ "- ' > 1 

it 



ir ' ^-W'Ml/J 9 '^"^""^' <28) 

This formula will be of crucial importance for computing hd(n, v, A), since it allows one to rewrite the 
integral in Eq. ([2|) as the following multiple integral over the manifold x §d-i x x §>d-i' 

hd(n, v, A) = 

—II / / / / (y^y d ~^' iK (y^r^- iv (y^r^- iX dya(^)a(^)a(^) 

(Zn) 2 io d _ x 

J lo J lo J la J <Ad 

(29) 

where we have used the measure (l26l) and the normalization ([23ll. 



2.3 Kallen-Lehmann-type representation for general two-point functions 

Consider again a general two-point function such that its reduced form u>(C) is analytic in the cut- 
plane II and uniformly bounded at infinity by a certain power |C| m °- It has been shown in [16) that 
for — 1 < too < there exists an integral representation of w(C) the following form: 

1 f°° ( d-1 \ 

w(C) = I I wo H 7. Mk I r (mo + d — 1 + in) T (— mo — in) x 

2(27r)2 J-oo \ 2 J 

x G(to + in) (C 2 - l)"^ 2 P"5L^ fC)d«5 (30) 



the function G(mo + in) is the boundary value of a function G(s) holomorphic in the half-plane 
Res > mo. G is obtained as a Laplace-type transform of the discontinuity Aro(() of u>(£) across the 
cut ] — 00, —1] (we will not need the explicit expression given in [IB], Eqs. Ill 10 and III 11). 

The results of [TB] can be extended to the case m = which is relevant for de Sitter quantum 

field theory, because, in that case, the Legendre functions involved in ([50)1 are all the free-field two- 
point functions of the principal series. We omit the details of the proof of formula ([50]) under the 
assumption that |itf(C)| is bounded by ICI - ^- Inserting the value too = — in Eq. ([30)) and taking 
the symmetry condition ([TT)) into account and puts , one obtains: 

poo Ip / d—l j A \ 1 2 d 2 i-OO 

w(C) = / Pi*) 1 \" ~ d n (C 2 " l)-^ 2 P:^ K (C)^ 2 = / p(«) Wk (C)^ 2 . (31) 
Jo 2(27r)2 2 j 

which is a genuine Kallen-Lehmann-type representation of w(Q) with weight 
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The computation the Kallen-Lehman weight p can also be tackled by invoking the generalized Mehler- 
Fock transformation theory [T7], as we do here. Eq. (f3"2")l takes the following concrete form: 

f 00 d — 2 

P (k) = 2 (2tt)^ sinh™ / w(Q PlPTJO (C 2 - 1)^ d(. (33) 



3 Decay of de Sitter unstable particles 

The study of particle disintegration in the de Sitter universe has been initiated in a pioneering paper 
by Nachtmann. 

Consider three independent neutral Klein-Gordon scalar fields (f>o, </>i, 4>2 with real mass parameters 
k, v, A respectively (i.e. the fields in the principal series) and an interaction term of the form 

7 g(x) £(x) dx, C(x) = : (f> (x)(f) 1 (x)(f)2(x) : 

where g is a smooth spacetime dependent " switching-on factor" which, in the end, should be made 
to tend to the constant 1. Self-interactions C(x) = : (f>(x) 3 : are a special case of this coupling. Let us 
consider the decay process 

0^1 + 2 (34) 
Let in particular \&o be a one-particle state of the form 

*o = / f{x)cl)o{x)fldx; 



the smooth test function f(x) contains the physical details about the preparation of the quantum 
state of the unstable particle whose disintegration we aim to study. The following general formula for 
the transition probability holds true [2J: 

7 2 27r|coth7rd f g(x) \F(x)\ 2 dx 

r(o ; i,2)= ' „J u y\ )\ y p v ,a(k). (35) 

Jf(x)W K (x, y)f(y)dxdy 

Here the convolution F(x) = J W K (x, y)f(y)dy is the "positive-frequency" solution of the KG 
equation with mass k associated with the test-function /; the denominator is the squared norm of 
^o- Note that the first factor in this formula does not depend on the the decay particles but only on 
the wavefunction of the incoming unstable particle. The infrared problem is contained in this factor 
and has to be overcome when letting the remaining g{x) tend to 1 (adiabatic limit). We will not treat 
this problem here and refer to [5]. 

The second factor is the relevant Kallen-Lehmann weight of the bubble diagram corresponding to 
two-point function of a composite field, obtained as the Wick product of the Klein-Gordon fields with 
mass parameters v and A: 

w(C) = M0wx(0- (36) 

This two-point function is well-defined and analytic in the cut-plane IT Moreover, for real values of 
A and v (i.e. for fields belonging to the principal series) it is bounded in II by |CI _d_1 an d therefore, 
a fortiori, by \C\ ~ . In particular, the Laplace-type transform G Ut \(s) of (|36j) is analytic in the 
half-plane {s E C; Res > — (d — 1)}. Thus, there exists a Kallen-Lehmann representation (j3T|) of 

w u (Qw\(C) = / Pu,\(n)w K (C,)dK 2 = I Kp u _ x (n)w K {()dK. (37) 
Jo Jr 
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The weight p, as given in Eq. (f32|) inherits analyticity properties from the aforementioned properties 
of G v x(s) of ([55)) and is itself holomorphic in the strip {seC: — (d — 1) < Re s < 0} and therefore 
it cannot vanish on any open interval of the line s = + £«■ 

This immediately implies that, in de Sitter spacetime, there is nothing such as the " subadditivity 
condition" of the Minkowski case: that property would require that p v ,\(n) should vanish if k < v + X 
and this would forbid the decay of a particle of mass k into a pair of particles of masses v and A when 
k < v + A. In contrast, in the de Sitter universe the disintegration of a particle of a given mass can 
give rise to two heavier particles if such a coupling enters in the interaction Lagrangian. 

In the following sections we will explicitly compute the Kallen-Lehmann weight by a mixture of geo- 
metrical insights and analytical techniques. By inserting the [16] of in the Mehler-Fock transform (|33[) 
of w v w\ it follows that p u ,\( K ) is proportional to the integral ([2]); more precisely, one obtains: 

P»,\(k) = '—^ ; \ 1 ) d 2 '—^ ^sinh(™) h d (k,u, A), (38) 

2(27r) i+ 2 



4 A special class of de Sitter three-point functions. 

The general properties of de Sitter two-point functions can, to some extent, be generalized to a class 
Cf of three-point functions W(xi, X2, X3) on de Sitter spacetime, such that 

1. W(xi , X2, X3) is a distribution on Xd x Xd x Xd which is decomposable as a sum of two boundary 
values of holomorphic functions W e (zi , z 2 , z 3 ), from the respective tubular domains T~ xT E x T + 
where e = + or — . 

2. Each function W e (zi, Z2, Z3) is invariant under the complex de Sitter group SOo(l, d)^ and 
therefore it coincides with a holomorphic function of the three complex invariants z, • Zj = 
dji hj = 1:2,3, (since zf = —1) in the image of the corresponding tubular domain. 

3. Each distribution W e (x\ , X2, x$) admits a Fourier- type transform on the one-sheeted hyperboloid 
defined in terms of the plane waves ([T9|) (see also [18] where the case d = 2 has been treated in 
detail) . 

We do not expect that the above properties hold for general interacting quantum field theories; in 
particular they do not apply to the general class introduced in [13] . They can be however useful in a 
perturbative context. 

The transform of W e (xi,X2,xs) is a distribution on (C + ) 3 x R 3 defined by 

W e {t 1 ,&,&;K,v,\)= [ (x 1 -Z 1 )+^- iK (x 2 -b)Z7 1 ~ i "(z3-&)Z^~ iX x 

xW e (xi, X2, X3) dxi dx2 dx3. (39) 

Conversely, each holomorphic function W e is recovered in its respective domain T~ x T e x T + by an 
inversion formula, which includes an appropriate weight- function ad on R + : 



W £ (z 1 , z 2 , z 3 ) = / ad{n)a d {v)<7d{X) dndvdX 

JR 3 , 



+ 

d— 1 



x / {z x - ^r^ +lK (z 2 ■ &)-^ +iv (z 3 ■ Z 3 )-^ +iX k, A)a(6)a(6)afe)- 



(40) 
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W e depends on (6,676) 0IUV through the (Lorentz) invariants 6 • 6, 6 1 6, 6 1 6, (since = 0). 
The homogeneity properties of (f3"9"|) w.r.t the variables 6 G C + imply that 

a d {K)a d (v)(j d (\) w £ (6,6,6;K 1 ^A) = p E (« )i /,A) (6 • 6) a3 (6-6) ai (6-6) a2 - (41) 

where we have introduced the parameters 

d — 1 ^ + A — k d — 1 A + k — f d — 1 k + i/ — A 
«i = 4 * ^ ' a2 = 4 1 2 ' 03 = 4 % 2 ' *■ ' 

The inversion formula (1401) can therefore be rewritten as follows: 



where 



W e {z ll z 2 ,z 3 ) = p e {K,v,\) w K , v .\{z l ,z 2 ,z 3 ) dndvd\ (43) 

JR+ 3 

w K>v , x { Zl ,z 2 ,zz) - / (^•6)"^ i+lK (^-6)"^ i+2!/ (^3-6)"^ i+4A x 

x (6 ■ 6) a3 (6 ■ 6) ai (6 • 6) a M6M6 W6)- (44) 

The three-point function w K! „ t \(zi, z 2l zs) manifestly satisfies the triplet of Klein- Gordon equations: 

[(KG) R ] Zl w K<Vj \(zt, z 2 , z 3 ) = [(KG) v ] Z2 w R ^\(z\,z 2 , z 3 ) = [{KG)x] Z3 w R ^\{z\,z 2 , z 3 ) = (45) 

in the (non-connected) complex open set (zi,z 2 ,z 3 ) e T ± x T ± x T ± where it is holomorphically 
defined via Eq. (|4"4"|) . This set contains in particular the relevant tubular domains x T e x T + 
in which the integral representation (|4"3"|) is meaningful. Formula (j4"3")l has the shape of generalized 
Kallen-Lehmann representation for all the three-point functions which belong to the class Cp on the 
basis of three-point functions satisfying the Klein-Gordon system (|451) . 



5 A star-triangle relation and a class of triangular invariants 
on the hypersphere 

Before the computation of hd(n, A) can be made possible we need to introduce two further ingredi- 
ents: a generalized star-triangle relation and a class of triangular invariants on the hypersphere. They 
both come out from the study of the following integral on the Lobachevski manifold y £ M. d (see Sec. 
[23: 

F 0l , a2 ,a 3 (6,6,6) - / (y6) Q2+Q3 (2/-6) Q3+ai (2/-6) ai+a2 %= (46) 

= c( ai ,a 2 ,a 3 )(6-6) a3 (6-6) ai (6-6) Q2 (47) 

the second equality follows again from Lorentz (i.e. de Sitter) invariance and from the homogeneity 
properties of (|46|) with respect to the variables 6's; the constant c(oi, a 2 , a 3 ) remains to be determined. 
This identity is a sort of "star-triangle relation" with one important difference w.r.t. what is usually 
called "star-triangle" : the center of the star is a point of M. d while the legs belong to the asymptotic 
cone C + i.e. the center and the legs of the star do not belong to the same manifold. By integrating 
both sides over the spherical basis 70 of the cone we get 

/(oi, 02,03) = / / / i? a 1 ,Q 2 ,a 3 (6 5 6 I 6)a(6)a(6)a(6) 

= c( ai ,a 2 ,a 3 ) I I I (6-6) Q3 (6-6) Ql (6-6) Q2 a(6M6M6) ( 4 «) 

•> 7o ^7o J~to 

= 0(01,02,03) X J(ai, 02,03) (49) 



Triangular invariants in de Sitter Universe 10 

The integral at the r.h.s. is proportional to the value of the three-point function iv K1/t x{zi, z 2 , 23) at 
the special complex event Z\ = z 2 = z 3 = iy £ T + , $/q = (1, . . . , 0): 

w K , v ,x{iyo^iy Q ) = e Ei!i ^ ±il e= 2ii ^ i [ (fc • 6)" 3 (6 • (6 • 6) aa a(6)«(6)«(6)- (50) 

•>7 3 

From Eqs (|29p. (|42l) and it follows that computing hd{n,v,\) is equivalent to integrating the 

star-triangle relation (|48|) w.r.t. the external legs: 

h d (K 7 v,X) = J 7(01,02,03) = 

(27Tj 2 

1 -7r(^ + ty + A) 3(d-l)7ri 

= - — -33 e 2 e w K , v ,x(iyo,iyo,iyo) X c(oi, 02,03). (51) 

(27T) 2 

Define 

7(01,03,03)= / &-6) M (6-6) Ol te»-6) 0B a&M&Wfc), (52) 

so that Eq. (|4"8"|) is rewritten as follows: 

7(01,02,03) = c(ai,a 2 ,a 3 ) J(ai,a 2 ,a 3 ). (53) 

The integral (|52p has a beautiful geometrical interpretation as a triangular invariant on the hyper- 
sphere [9] Sd-i. Consider indeed the squared distance Arr? fe between two points n^ and n& belonging 
to §d_i. The Lorentzian scalar product of two points £j = (l,nj) and = (l,nfc) belonging to the 
spherical cycle 70 of the forward lightcone C + is proportional to the squared distance An? fc : 

An 2 ik = (m - n fc ) 2 = 2 - 2 m ■ n k = 2& • &. (54) 



Given three points ni, n 2 and n 3 of §d-i and three complex numbers Oi, a 2 and 03, we construct the 
rotation invariant quantity 

J = .7(01,02,03) = ((AnL) a3 (Anl 3 ) ai (An^f 2 ) (55) 

where (/) denotes the average on §<j_i; for three points 

L 3 dnidn 2 dn 3 /(ni,n 2 ,n 3 ) 

(/> = -^7 , , , • (56) 

J S 3 i dni, dn 2 , dn 3 

It follows that 

J(oi,o 2 ,a 3 ) = ^ / (a-6) a3 (6-C3) ai (C3-Ci) a2 ^7D(6)^7o(6)rfM7o(6) = ^J(oi, a 2 , a 3 ) (57) 

and therefore, in view of (|53|) : 

7(01,02,03) = 2 _a ^ 0(01,02,03) J(ai, a 2 , a 3 ) (58) 



In the following sections we will describe the details of the concrete evaluation of hd(K, v, A). The 
method that we present below is based on the previous relations (|4"5|) and (fijTj) applies to all values of 
the spacetime dimension d (d > 2). 

There exists a remarkable symmetry relation between of c and J that follows from symmetry rela- 
tion 

hd(K, v, A) = hd(—K, —v, -A) (59) 
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which in turn is a consequence for all real values of k, v, A, of the basic symmetry property (|17p of the 
Legendre functions first-kind. If we introduce the corresponding transformation 

, d- 1 

a i -> a i = ~ a J — > 3 = 2 ' 3 ' ( 60 ) 

and take Eqs. (|42p and (|51|) into account, we obtain that for every triplet (01,02,03) such that 
Reoi = Re a- = — =ji, the following equality is valid 

/ (01,03,03) = /(oi, a 2 ,a 3 ) (61) 

or either (in view of Eq (|55|) ): 

2 _ a J(ai,02,a 3 ) = ^ JK.fl;,^) , g2 ) 
c(ai,a 2 ,a 3 ) c(ai,a 2 ,a 3 ) 



which are both equivalent to ([5 

This striking duality between two integrals over different manifolds is surprising at first sight. It will 
be made clear that the factorization corresponds precisely to a splitting of the expression © of 
hd(K, v, A) into two parts which are symmetric under that parity transformation. 



In the following section we will evaluate the functions c and J. The study of the integrals (146|) and 
(|52|) will lead us to define and compute the functions c and J in appropriate domains of the complex 
space C 3 of the variables (01,02,03). hd(K, v, A) is then obtained by taking the restriction of the 
holomorphic function / = c x J to the linear real submanifold Ld of C 3 defined by the equations (|4"2")) , 
namely L d = {(01,02,03) € C 3 ; Reoj ; = -^=±; j = 1,2,3}. 

6 Computing c(a\, 0,2, 03); more on the star-triangle relation 

For computing the integral (|46[) and obtaining (|47|) with the complete expression of the function 
c(oi, a 2 , 03), we shall consider the following double Mellin transform: 



OO /-00 



JO 



T(-b) 



this relation is valid for Re(s) > 0, Re(i) > 0, Re[& + s + t] < 0. Mellin's inversion theorem then 
provides the following expansion: 

, b _ i r +MO , - r' +icc ^^ t ns)T{t)T{-b-s-t) 



(l + v + z T = -— dsv-s dtz- — (64) 

J ioo J 7 —zoo L V u ) 

where the integration paths lie in the strips allowed by the previous inequalities. The formula can be 
rendered symmetric by homogeneity: 

1 r +i °°^ r ,+ico A „^t..-..-trWr(t)r(-t -s-t) 



(u + v + z) b = -— ds dtu b+s+t v- s z~ t v ; v ' v '-. (65) 

J~y — ioo J ~{ f — ioo \ V) 

Consider now the expression (y ■ 2) 2a where 

5 = ^+6+6, a = ai+a 2 + a 3 ; (66) 
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all the £j's are lightlike and therefore S is either timelike or lightlike (in the latter case £1,^2 and £3 
lie on the same generatrix of the cone): 

S 2 = 2£i-£2 + 26-£ 3 + 2£ 3 -£i>0. ( 67 ) 
By application of (f6"5)) we get that 

(r^ = -^ r°° ds r' +i °° dt r(s)m rtt~ 5 ~ V ^) 2a+s+t (y ■ br B (y • &r* (68) 

with 7 > 0, 7' > 0, 2 Re(a) + 7 + 7' < 0. 

Let us integrate the two members of (155]) over the Lobatchveski manifold H^. The Lorentz invariance 
of the l.h.s. implies that the integral can be parametrized in terms of the hyperbolic angle A between 
y and 5: 

(y ■ Zf"dy = -^(S 2 )" H (cosh af« (sinh a^da = ° " fi (5 2 )° i (69) 

this result holds provided Re(a) < — (d — l)/2. Thus 

2 a vrir(i-a-|) _ 1 2 a 7rir(i-a-f) 

pzjr^ + --4^2 r (i_ a ) r( _ o) x 

ds / *(a-6) Q+s+t (6-6)~ s (6-6)^rwr( s )r(-a- s -i) = 



8—ioo J 8' — ioo 



1 



7+zoo />7'+'too 

ds / dt - 



J^i-ioo J 7' -ioo r( 2a) 



(70) 



In the second step we have applied once more Eq. (|65|) with S > 0, (5' > 0, Re(a) + <5 + 5' < 0. Since 
Eq. ([7D]) is valid for any choice of the vectors £1 , £2 and £3 we can multiply £2 and £3 by two complex 
numbers a and (3 and obtain that 

7+ioo /'7 / +'too /> 

ds«- s / dt/r<r(*)r(t)r(-2a-»-t) / d 2/ ( 2/ .£ 1 ) 2a+s+t (2/-6)^(y-6)"* = 

-.2- a -^r(--a--)x 

dua a+u / d« /3- u ~ v T(u)T(v)r(-a - v - u)fa • 6) a+l ' + "(6 • 6)^(6 ■ 6)"" 

(5— -ioo J 8' —ioo 

(71) 

By changing the variables in the second integral as follows: u = —a — s, t> = a + s + t. the r.h.s. 
becomes 

- Rc(a) — (5+ioo +ico 

ds a~ s / dt /? _t r(-a - s)r(a + s + f)r(-a - t) x 

R.e(a) — 8 — zoo J 8-\-8' — ioo 

x(£l-6) 0+t (6^3)-°- S -*te-£x) a+S . 

(72) 

By Mellin's inversion theorem we can now identify the integrands i.e. 

'*/(y-£x) 2a+s+ V6rV6r' = 
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2-«-VgT (|-a-j) T(-a - sjTja + s + t)T(-a t) _ 
rmtn _ 2a - s -t) & ■ & ' 6) & ' a) 

(73) 



Finally, by setting — s = 03 + ai and — < = ai + 02 the proof of the star-triangle relation (J47J) is 
completed with 

„ 1 d-i_ ( d — 1\ r(— ai )r(— a9)r(— a.s) 

C (a 1 ,a 2 ,a 3 ) = 2- a -V^r - a - __ 1 * 1 V \ r( 3 ' r (74) 

\ 2/1 (-02 - 03 )1 (-03 - 01 )1 (-ai - a 2 ) 



6.1 Conical limit 

It is interesting to remark that, under certain conditions, a true star-triangle relation is obtained by 
integrating over a cycle of the lightcone. Consider indeed the three-point function 

G ai , Q2 , Q3 (6,6,6) = / (6 6) a2+a3 (6 6) a3+ai (6 6) Q1+a M£), (75) 

where the integration is performed on the parabolic basis C + fl {£° + ^ d = 1} of the future cone; this 
corresponds to setting A = 1 in the parametrization 




1(1 + v 2 ) 



Z(\v) = \ € = W 0<A<oo, !)£IR^, (76) 

= fa-'? 2 ) 

and w.r.t. the Lebesgue measure d/i(^) = drj. Concretely since 



2£(A,7 7 ).£'(l,r/) = A( ?? -7/) 2 (77) 

>02+»3 \ «3+»l \ "l+"2 P 

G aua2 , a3 (6,6,6) = ^ |^ y [(^-'7i) 2 ] a2+a3 [(^~'72) 2 ] a3+ai [(^-'73) 2 ] ai+a2 ^. 

(78) 

The result is of course invariant under Euclidean transformations in (d — 1) dimensions. There is 
however a special case to be considered: when 2a = (1 — d) the integrand is a closed (d — l)-form 
on the future cone. By applying Stokes' theorem one sees that the integral does not depend on the 
choice of the integration manifold and the result is fully 6*00(17 (i)-invariant; this invariance can be 
now interpreted as Euclidean conformal invariance. As before, by exploiting the SOo(l, d)-invariance 
of the integral and the homogeneity properties of the integrand one has that 

K • £i) a2+a3 (£ • 6) Q3+Q1 (? • 6) ai+a3 <M0 = d(oi,aa,03)(6 • 6) Q3 (6 • 6) ai (6 • 6) Q2 - (79) 

M 

To determine the constant we need to compute the integral J (£ -y) 1 a(£) where again y = 6+6+6- 
This is most easily done using the spherical basis of the cone 70 with respect to the measure rf/x 7o (£) 
(see Eq. (|23"]) ). Calculating the previous integral in these coordinates is immediate: 

it'V) 1 - d <*(£) = o>d(yv)- s ? (80) 

Now we can follow the same steps as before and obtain the limiting conformal invariant star-triangle 
relation (75]) with 

d(a 1 ,a 2 ,a 3 ) = (2tt) s> — — — (81) 

1 (-02 - 03)! (-0,3 - ai)l {—ai - a 2 ) 
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7 Computing J(ai, 02,03): direct and indirect methods. 



7.1 Probabilistic interpretation of J 

To evaluate J we make use of the standard parametrization of a point n G S<j— 1 in terms of (d — 1) 
angles 0i, parametrization that we spell here for the reader's convenience: 



n = < 



cos 6\ 

sin 81 cos 6*2 

sin 61 . . . sin 64-2 cos 0^-1 
sin 0i . . . sin 64-2 sin 9d-i 

so that 



0< 6>l <7T, ...,O<0 d _2 <tt , s 

< d _i < 2ir [SZ) 



dn = Y[(sm6 j ) d - 1 - j dO v (83) 



The rotation invariance of (|55[) implies that it is possible to perform the integrations by fixing 
" point, say ni = (1,0,.. .,0) and specializing the second point n2 = (cos cf)\, sin <pi, 0, . . . , 0) so 



one 
that 



n.1 • n.2 = cos <fii, ni • n3 = cos 0i , n2 • 113 = cos <f>i cos 9% + sin <pi sin Q\ cos 02 ■ (84) 
It follows that 

J = *^zpk± [* (wxfa)*-*d4i r {sm6 x ) d - 2 d0i (sin e 2 ) d - 3 d6 2 x 
Jo Jo Jo 

1 — cos 0i \ 03 / 1 — cos 0i \ ai / 1 — cos (j>i cos 0i — sin 0i sin 0i cos 2 



■ (85) 

Before proceeding with the evaluation of J, let us change the integration variables at the r.h.s. of 
(f85|) and replace the angles used there with the distances 



ri=An 2 3, r 2 = Ani 3 , r 3 = Ani 2 , < r,- < 1. (86) 
The Jacobian of the transformation is readily computed: 

riT2T3 dri dr2 dr 3 — (sin^i) 2 (sin0i) 2 sin 62 d(f>i dd\ dQ 2 (87) 

Since 

(sin 0i sin 0i sin0 2 ) 2 = 2(r 2 r 2 2 + r\r\ + r\r\) - {rf + r\ + rj) - r\r\r\ 

it follows that 

J(ai,a 2 ,a 3 )= / p(n, r 2 , r 3 ) r\ ai r 2 " 2 r 2a3 dridr 2 dr 3 . 



D 



p(ri,r 2 ,r 3 ) = n r 2 r 3 [2(rl4 + r\r\ + r 2 r 2 ) - (rf + r\ + rj) r^r 2 ]^ . (90) 

(1} . 



where 

4 Q Wd-i^d- 
Below it will be proven that 

J p{r 1 ,r 2 ,r 3 )dr 1 dr 2 dr 3 = 1 (91) 

J(ai, 02,03) are therefore the moments of the probability density of three random points on S^-i 
constituting a triangle whose sides have the sizes n, T2 and r 3 . 
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By introducing a variable z as follows: 

1 — cos 4>i cos 9i 

I COS 0i — COS 0i I 



vV-i 



sin 0i sin U\ 
I cos 0i — cos 6>i I 



one can identify a well-known integral representation of a Legendre functions of the first kind entering 
at at the r.h.s. (see e.g. [T] Eq. 3.7.6): 



J = 4 q (2tt 



d-l LUd-l 



d JO 



x 2 



d-3 



COS (pi — COS tfi 



d cos 0i / <i cos 9i 
Jo 

d— 3+a2 



1 — COS 01 \ ' / 1 — COS 



(92) 



Evaluating this integral involves several nontrivial steps. In the first one we change to the variables 
x = (1 — cos0i)/2 and y = (1 — cos6>i)/2 and restrict the domain of integratation to the region 
< y < x < 1: 



d JO JO 



dw(z - w)- a2 -\w 2 - l)^+ aa ; 

(93) 

in this Equation we have made use of the following remarkable integral representation of P„{z) that 
will be established in Appendix [A"l 



(z>-1)^P- > + Jz) = 2- 



dw 



(z - w)- a2 - 1 (w 2 - 1) — + aa 



i r(-a 2 ) r(^i + a 2 



(94) 



In the second change of variables we replace y by [ix: 



J i = 



4°+W^(2,)^Wi ai(1 _ ^-s.,, 



-+a 2 



g?u>- 



"7 \ 1 1 : (95) 



2 can be expressed in terms of x and /x as follows 



z — z^ -\- (1 z^)x^ z^ 



1 + jH 



(96) 



Next, we use z as integration variable in place of x (at constant [i) and we get that 



dxx a+d - 2 



dw {z-w)- a2 -\w z - 1) 



Zfj, - 1 



i V I+ ' i " 1 r(a + rf-i)r(- a2 ) 



r (a + d — 1 — a 2 ) 7i 
By setting u = (z^ — w)/(z^ — 1) we obtain the final expression for Ji 

A a 2^(2iT)^ 1 u Jd -2 r(o + d-l) 



dw (Z M - u ,)«i+«3+d-2( u; 2 _ !) — +a 2 _ ( 97 ) 



r(^ + a 2 )^_i r(ai + a 3 + d-l) 



dm^ +a3 (l-u)^ +a2 / v^+^(l-v)^ +a2 dv 



(98) 
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The integral J 2 over the domain 0<x<y<lis obtained from this expression by interchanging a\ 
and 03. Finally 

T , 7 2^+ 2o (2 7 r)^ i o; d _ 1 r(o + d-l) 



(99) 



In the end 



r(^ + a 2 )^ r( ai + a 3 + d-l) 

2 ^+2a (27r) ^i CJd _ i r( a + d - i)r (jg + ax) r(jg + a 2 )r + o 3 ) 

r (d - 1 + a 2 + a 3 ) T (d - 1 + a 3 + ai) T (d - 1 + ai + a 2 ) ' 



h d = /n \ M 2 a ujjc(ai,a2,a 3 )J(ai,a 2 ,a 3 ) 

(27T) 2 

2 ^r(- Q - } r( a + d-i) JL r + a,) r(- aj ) 

(27 r)3/2r(^i) llrfd-i + o^-ajr^-o) 1 j 

and the formula ([2]) is proven for d > 1 integer. Estimates shown in Appendix JB| show that 
interpolation is permitted and the formula is vaild for complex values of the dimension d. The 
analyticity in k of the expression at the r.h.s. also directly confirms our previous statement the 
analyticity of Pv,\{k)- 

7.3 The J/c— duality and an indirect evaluation of J(ai,a 2 ,a 3 ) 



There is a possibility to avoid the computation based on the relation (|62|) and on the fact the inte- 
gral ([57]) defines ,7(01,02,03) as a holomorphic function of (01,02,03) in the domain {(01,02,03) £ 
C 3 ; Re(ai + a 2 + o 3 ) > -(d - 1)}. 

It then follows from the latter property and from the analyticity properties of the T— functions involved 
in the expression J73|) of the function c(a[, a' 2 , a 3 ) (with = — Oj — that the l.h.s. of (p)2"|) can 
be analytically continued as a holomorphic function £(01,02,03) of (01,02,03) in the tube T + = 
{(oi, 02, 03) £ C 3 ; Reai > — i = 1, 2, 3}, while the r.h.s. of (f6"2"|) can be analytically continued as 
a holomorphic function r(ai, 02,03) in the tube T_ = {(01,02,03) £ C 3 ;Reai < ~^pS * = 1,2,3}. 
Note that the two tubes T + and T_ , which are symmetric of each other with respect to the submanifold 
Ld, have a nonempty intersection which is a tube containing Ld- Since both sides of (f6"2"]) coincide on 
Ld, it follows that they both represent the same analytic function -E(ai, 02, 03) in the complex domain 
T + n T_ , which therefore also admits an analytic continuation in T + U T_ . But in view of the tube 
theorem, the holomorphy envelope of T + U T_ is the full space C 3 , and thereby _E(ai,02,a 3 ) is an 
entire function on C 3 . 

Moreover, one can show that the functions I and r are uniformly bounded by a constant C in respective 
tubes T' + and T'_ , which are conical open neighborhoods of the respective "diagonal tubes" <5 + = 
{(oi, 02,03) £ C 3 ;Reai = a > — 2^1; i = 1,2,3} and S- = {(01,02,03) £ C 3 ;Reai — a < 
— i = 1,2,3}. Since the convex hull of T' + U T'_ is equal to C 3 , it then follows that the entire 
function E{ai, a 2 , a 3 ) is also uniformly bounded by C in the whole space C 3 and is therefore a constant 
function. This entails that the l.h.s. (or r.h.s.) of (p)2"|) is a constant 7^ (only depending on d) and 
that the integral (|5T[) is directly obtained without additional computation by the formula 



7(01,02,03) =7^ 2 a c(ai,a 2 ,a 3 ) (or J(a[, o 2 , a 3 ) = j d 2 a c(ai, 02,03)) 



(101) 
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Since J(0, 0, 0) = 2 a , the constant jd is easily computed by choosing aj = in (jlOip . which yields (in 
view of (|71)) ): 

3d- 5 d-i [T(d- 1)1 2 . , 

u [r(^i)] 3 v j 

and therefore (in view of (158|)): 

w 3 2 2d ~5 

h d (K,v,\) = d d 7 d £(0^,02,03) x c(a 1 ,a,2,a 3 ) = , cfa^, a 2 , a 3 ) x c(ai,o 2 ,o 3 ) (103) 

(2tt) 2 7r a l ( 5 ) 



7.4 A corollary 

A formula by Hsu quoted in Szego's book ([TS], page 390) gives a weighted integral of a product of 
three Gegenbauer polynomials, all having the same upper index k — That formula is a sort of 

extension to SO(d) of a well known relation established for Legendre polynomials (upper index k = 1) 
and expressing the square of a Clebsh-Gordan coefficient of the rotation group 50(3). 

The Kallen-Lehmann formula given in Eqs. (|3T[) and (|38[) is an identity between two functions 
holomorphic in the variables £ and d. It is possible to perform an analytic continuation of (|37[) to real 
values of the variable £ belonging to the interval ] — 1, 1] for every d. 

By studying the behavior at infinity of the integrand one can show that the integral converges uni- 
formly w.r.t. £ in that interval (with the exclusion of the singular point £ = —1). Then, it is tempting 
to compute the r.h.s. as an infinite discrete series of residues (there is no difference in closing the 
contour to the right or to the left). But the so obtained series converges only for k < 1 i.e. d < 3. 
The residue of the pole in in = k + n (with k < 0) of w K is the Gegenbauer polynomial of degree n. 
By specializing the l.h.s. of ([3T]) to iv = k + n\ and iX = k + 122 the series at the r.h.s. becomes a 
finite sum in the range \m — n 2 | < n < \n\ + n.2 1 - The aforementioned formula by Hsu and Szego is 
thus a corollary of our result ([2]) via the analyticity of the Kallen-Lehmann formula (|37| . 



8 Conclusions and outlook 

The integral {l} that we have studied in this paper gives an exact evaluation of the "phase space" 
coefficient in the rate (|35[) of the decay process 1 — » 2 in a de Sitter universe. The fact that p v ,\(n) 
is an holomorphic function in their arguments means that in the de Sitter universe all such decays 
are possible and there is no analogue of the mass subadditivity condition of flat spacetime. This 
phenomenon should not be ascribed to the thermal features of the de Sitter universe [T51 [TTJ [TJJ 

na. 

Indeed the thermal features manifest themselves only when the field algebra is restricted to a wedge- 
shaped region bordered by bifurcate Killing horizons 20J while the amplitude (|33|) is integrated overall 
the de Sitter universe. Second, a similar computation performed in flat space thermal field theory 
does not exhibit the same phenomenon. 

There several possible other investigation that are rendered possible by the results of the present 
paper. One is the exploration of the star-triangle relation that we have found and the associated 
Yang-Baxter equation. We will examine that question elsewhere. 
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A Appendix. Another integral representation for the asso- 
ciate Legendre functions of the first kind 

In this appendix we show how to construct the integral representation of the Legendre function of the 
first kind that enters crucially in our computation of the triangular invariant ,7(0,1,02,(23). 

The starting point is the usual expression for the associated Legendre functions as derivatives of 
Legendre polynomials commonly encountered in books of quantum mechanics: 

lrn 1 jm+I 

Pr{x) = (x* - = W y - l)-^Tr(- 2 - 1)' (104) 

valid for x > 1. The idea is to use the Riemann-Liouville fractional integration operator 

(D v f)(x) = — ^ J\x - y)-^ 1 f(y)dy (105) 

to replace the standard derivative of integer order in (|104[) . The conjectured interpolation formula for 
the Legendre function of first kind is then given by 

™ - = 2>r(A ( +iffH-A) f < z - »>-"*" V - < 106 » 

To directly evaluate the RHS of (| 1 06[) we set z — 1 = 2£ with — n < Arg£ < ir and y — 1 = 2A£ with 
< A < 1: 

= JrP-T^^T, / 1 (l-A)-^ A - 1 A A (l+CA) A dA 



r(A + i)r(- M -A) J 

L(A + n+ 1) 



= (z-iyHz + l)^? 



n=0 



r(-/i + n + 1)L(A - n + l)T(n + 1) 



( g _l)-f (x + 1 )f ~ r(A + n+l)r(-A + n) 

tt ^ ^ r(-M+n+l)r(n + l) 

n— J 

- {z ~ ^ f _ ( ;+ 1)f ,Fx (-A, A + 1; 1 - m; ^) = ni*). dor) 

The last identification follows from [1] Eq. 3.2(14). The proof is valid in the domain of convergence 
of the hypergeometric series but the identification has of course a larger domain of applicability that 
is the cut-plane C \ {(-co, 1]} where the integral at the RHS of (|106|1 converges. A second related 
way of evaluating () 106(1 consists in using Mellin's convolutions and writing 

(x - J,)""-*- V - l t d y = x^- 1 r f (-) 9(y)~ = x-^-'F * G(x) (108) 

1 Jo w J y 



- \ -n-X-l 



where 

f(x) = (l - , g(x) = x(x 2 - 1)\. (109) 

The Mellin transforms F(s) and G(s) of these functions are readily computed and their product is 
the Mellin transform of the convolution (|108[) 

F(s)G(s) = 2 A +^r(-/x-A)r(A + i) , r H)r(-f-^-A) 

£ 4. ~ M~ \ p I _ a 1 — ^— A+j 

2 ' 2 y ^ 2 ' 2 
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Inversion of the Mellin transfom can be obtained by applying Slater's theorem [21] or either directly 
by integrating a la Mellin-Barnes; after a few steps that we do not reproduce we get again that 



F * G(x) = — / F{s)G{s)x~ s ds = 2 X T{\ + l)r(-/x - \){x 2 - l)-%P£(x) (111) 
and the representation is proven to hold true. 

There exist of course many other integral representations of Legendre functions. The advantage of 
this specific one is that it can be used easily for partial integrations and is very natural generalization 
of the classic expression (|104|) . Strangely enough, we were not able to find this explicitly written in 
the form given in Eq. (|106p in the literature available to us. 



By using the integral representation (|106p for real u > 1 we get the following bound 

n 1 , _ . Rc d i 1 , _ . Rc d i 

23 (u + 1) 4- + 3(w _ I)— 1 



< 



|r(!fci-ii/)r(i + ii/)| 



log (u + \fl. 



(112) 



so that 



h d (n, v, A) 



< 



2iC d (K, v, A) + l)- Red + 2 (tt - 1)^ [log (u + du 

|r(i + iAt)r(i + ^)r(i + iA)| 
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The change of variable u — cosh</> allows to rewrite the integral at the RHS of Eq. (|113[) as fol- 
lows: 
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(115) 



The asymptotic behaviour at large d of the integral at the RHS has been estimated by the the 
steepest descent approximation; for large Re d the maximum of the integrand is located near to 
4> c = 2 log (l + \/2) . The behaviour of Cd{n, v-, A) can be obtained from Stirling's formula (PQ , p. 47) 
that tells us that it goes like a constant. It follows that 



< const - 



2~3Re d/2 



Vied 



(116) 



This asymptotic behavior will be promoted to a bound in the following section. 



B Appendix. Extension of the main formula to complex d 

We have shown that the function h d (K, v, A) defined in Eq. ([1]) coincides with 

, x] = (4,)fr(^i) U + 2 ) 

9d[K ' ) [IL= ± i r + <«)] [iu ±1 r + ie> v )] [n e »= ±1 r + ie"A)] 
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(117) 

for all integer d > 2 and real n, v, and A. In this appendix, we will show, by using Carlson's Theorem 
|22j . that these two functions coincide wherever they are both defined, i.e. where the integral in Q16p 
converges. It is obviously sufficient to prove this for real values of k, u, and A satisfying \k\ < B, 
\v\ < B, |A| < B for some arbitrary B > 0. 

Theorem B.l (Carlson) Let f be a function holomorphic in the right half-plane {z S C : Rez > 
0}, and satisfying 

l/OOl < Ae m for all z with Rez > 0, (118) 
where A > and < k < n. If f(z) = for z = 1, 2, . . then f = 0. 
Recall the formula (pQ, 3.7(1) p. 155) 

P£(z) = W, J I TT / (^ + cht)^- 1 (sh<) 2Q + 1 *, (119) 



r(- M - a )i> + 1) y 

valid for z G C \ (—00, 1] and Re(— /i) > Re a > —1. We apply this to the special case Re a = —1/2, 
which satisfies the above condition when — Re/i > —1/2. We shall suppose — Re/j, > 0. In this case, 
for z > 1, 



|2 a r(-^-a)r(a + l)||P£(z)| < (z 2 -l)- Rc ^ 2 / (z + chtf ^ 1 / 2 dt 

Jo 
1*00 

< {z 2 -l)- Rc ^ 2 (z + l^-^+^chi)-^ 

Jo 

poo 

< (z 2 - i)-B»ii/2^ + 1 jltoM-i/a+s / 2 £ e~ £t dt 

Jo 

< 2 e e- 1 (z 2 -l)- Rc ^ 2 (z + l) RcM ~ 1/2+e (120) 



This holds for all real e such that < e < — Rc/i + 1/2, in particular for < e < 1/2. Note that 
setting e = (2 + logz) -1 gives a more precise bound, but (|120[) will suffice for our present purpose. 
Keeping k, v, and A real and setting /i = — [d — 2)/2, with Red > 1, this gives 

2 |+3e £ -3 

M ^' A)I " |r(^-z K )r(i + iK )r(^i-z,)r(i + w )r(^i-zA)r(i + zA)| x 

(z 2 - l) Rc ^ (z + l)- 3Rc ^ +3e dz (121) 



Bounding (z — 1) by (z + 1) , the last integral is bounded by 

2 (Re^i-3e) V ; 

provided 3e < Re(d - l)/2. We choose e.g. e < 1/12. Thus, for Red > 2, 
|/i d (t, j/, A) I < 

2§+6s-Rc(ci-l)/2 e -3 

(Re 3=1 - 3 e ) |r - m) r (| + i«) T - iv) F(f + if) f - iA) f(i + iA) | 



Recall Stirling's formula (PQ, p. 47): 

T(z) = (27r )* e -*+C*-*)>«** (1 + ^ + £_! + 0(z -3)^ 



(123) 



(124) 
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valid for z£C \ R . In applying (] 1 24[) . we will ignore the last bracket, at the cost of introducing an 
unknown multiplicative constant and supposing \z\ > a > 0. Thus there is a constant K such that, 
for z = x + iy, Re z > and \z\ > a, 

\r(z)\<K\zr^ 2 exp(-x + x\og\z\-^\y\), ^(z)]- 1 < K\z\^ 2 exp(x - x\og \z\ + ~\y\) . (125) 

If |«| < B, |r(| + * K )| _1 is bounded by a constant (depending on B). Assume Red > 2 and 
\d — 1| > 2e + 2£>. Applying the second inequality in (|125| with z = (d — l)/2 — zk we find that, for 
some if' > 0, 

|r(^-4r(*-H«)| " KV - 1,1/2 exp(7rM - 1|/4) • (126) 

Therefore, for Red > 2, \d - 1| > 2e + 25, |k| < £?, \v\ < B, and |A| < B, we find that, for some 
constant A\ > (depending on £>), 

v,X)\< A 1 \d - i\m e Md-W . (127) 



With the same assumptions on k, v, and A, we now wish to apply (|124[) in formula (jTTTJ) , taking 
advantage of the cancellation of growth which occurs between numerator and denominator. Let 
z = x + iy = (d — l)/2. The function g d has the form 



g d (K,v, A) = 22(4tt) = 



nr(i 

J'=l 



]7r(z + zt- fc ) 1 



Lfc=o 



(128) 



where u.,-, Wfc are real, |tt 3 -| < 25, < i3, (vq = 0). We suppose x > 1, |z| > 4B + e, \y\ > 4B. There 
is a constant > such that 



\g d (K, u, A) 
exp 



< A 



2\Z\ 



f(iogN +iogia +1^1/21*11 -1)- 



i=i 



2B 



+ 



u 

^ z(- log \z\ - log |1 + iv k /\z\\ + 1) + + B) 



.fe=0 

With our assumptions, log|± + zMj|/2|z| | < and — log 1 1 + Wfc/|z|| < 0, hence 

\g d (K,v,\)\ < A 2 \z\~ 1 exp 
and finally, for some A 3 > 0, 



3a;(log \z\ - 1) + + 15ttB/2 



\g d {K,v,X)\ < i4 3 |d-l|- 1 exp 



37r|d- 1| 



(129) 



(130) 



(131) 



Therefore the function d <— > h d {n, v, A) — g d (n, v, A) vanishes by Carlson's theorem. 



C Appendix. Analytic continuation of the Kallen-Lehmann 
weight 



In this appendix, we use analytic continuation to obtain the Kallcn-Lchmann weight for a two- 
point function which is the product of the two-point functions of two free fields belonging to the 
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complementary series (a special case was done in [3]). We fix R = 1, and d > 2 is an integer. We 
have obtained the representation 

Wu{()w\(() = / kp„,\(k)w k {C) dK , (132) 

with 

1 K sh7TK 

Kf>uM = 2VTO^r (^i) r + i«) r - i«) x 

x jj r + i^±i^+i^ \ . (133) 

e, e', e"=±l ^ " / 

Eq. (| 1 32[) has been proved for real values of f and A. We know however that the l.h.s. of this 
equation extends to a meromorphic function of v and A with simple poles at %v — ±(2=i — n) and 
iX = ±(=7p — n.), where n > is an integer. As i/ and A become complex, (|132l) will remain valid 
as long as the r.h.s. can be continued. The integrand in the r.h.s. is a meromorphic function of k, 
v and A. We denote /i = ==I, In the variable K, the poles arise from the last product of Gamma 
functions. When their arguments arc close to —n, where n is any integer > 0, these T functions 
behave as follows: 

r(n+~(K + e'p + e"\)) 2 }~ l) l., —j^ (134) 

\ 2 / nh\K — i(2(i + 2n + te v + le A)) 

r( M +~(- K + eWeV) ~ — ~ 2 [~ 1)n . , — — (135) 

V 2 V ') nh(K + i{2p + 2n + ie'v + ie"X)) 

Thus the poles listed in (|134[) are the opposites of those in (|135D . 

When v and A are real, and more generally when | Im^| + | Im A| < 2/i, these poles do not touch the 
real axis, so that (|132[) holds, by analytic continuation, for all such values. 

For definiteness we first assume < j3 = Im A < a = Im v. We also temporarily assume that Re v =^ 0, 
Re A 7^ 0, and Re v ± Re A 7^ 0. The poles are at 

en = -e' Re v - e" Re A + i(2fi + 2n - e'a - e"(3), (136) 

They can be grouped as follows: 

1. On the line — Rc v — Re A + iR : 

(e=l, e' = e" = l) : k = - Re v - Re A + i(2p + 2n - a - (3) ; (137) 

(e= -1, e' = e" = -1) : n = - Rc v - Re A + %{-2\i - 2n - a - f3) ; (138) 

On this line, all poles move down as a + (3 increases, The mutual distances of these poles remain 
constant. 



2. On the line Re v + Re A + iR : poles opposite to the preceding. 

3. On the line - Re v + Re A + iR : 

(e = l, e =-e" = l) : n = - Re v + Re A + i{2 f i + 2n - a + 0) ; (139) 
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(e = -1, e' = -e" = -1) : k = - Rev + Re A + i(-2/i - 2n - a + /3) ; (140) 

On this line, all poles move down as a — (3 increases, The mutual distances of these poles remain 
constant. 

4. On the line Re v — Re A + iR : poles opposite to the preceding. 

When a + (3 increases and reaches 2/i, the pole in (|13T|) with n = reaches the real axis, as well 
as its opposite. Beyond this, eq. (|132p aquires two discrete contributions provided by the residues 
of these poles; these can easily be seen to be equal. If a + (3 increases further, poles in (|137[) with 
higher values of n cross the real axis, contributing new discrete terms to eq. (|132l) . In the range 
2/j + 2N < a + f3 < 2/j + 2N + 2, (with integer N > 0), all the poles in (fl3T)) with < n < N have 
crossed (as well as their opposites). 

Similarly if a — f3 increases so that 2/i + 2N < a - (3 < 2[i + 2N + 2, (with integer N > 0), all the 
poles in (|139p with < n < N have crossed the real axis (as well as their opposites). 

Therefore, if 2fx + 2N < a + (3 < 2fi + 2N + 2 and 2fi + 2M < a - /3 < 2/x + 2M + 2, where M and 
N are non-negative integers, and with our assumptions regarding Re v and Re A, 

w u (z, z')w\(z, z') = / np v .\{n) w K (z, z') dn 
Jr. 

N M 
+ y^ j A n (u, X) Wi( 2 »+2n+iv+i\) + X! B n( U > w t(2^+2n+iv-i\) ■ (141) 



Here 



n=0 n=0 



M", A) = — d+1 



n!2 2 ?rT-i? d - 2 r (^=i) T(-2fi -2n-iv- «A)r(2/x + 2n + iv + iX) 
1 



x 



and 



r(-2n -iv - i\)Y{Ap, + 2n + iv + iX) 
]?(— n — iv) r(— n — iX) T{— n — iv — iX) x 

r(2^i + n + iv + iX) T(2fi + n + iX) T(2fi + n + iv) T(2p + n) , (142) 



B n (v, A) = A n (v, -A). (143) 



In the expression for A n {y, A), we note that the poles in T(—n — iv — iX) are cancelled by those of 
r(— 2n — iv — iX) in the denominator. The poles of T(2p + n + iv + iX) are cancelled by those of 
T(2/i + 2n + iv + iX) in the denominator. Hence the possible poles of A n (v, X) come from the poles 
of: 

r(— n— iv) : —n — iv = —m Rc^ = 0, a = n — m ; 

r(— ri — iX) : —n — iX = —m RcA = 0, [3 = n — m ; (iaa\ 

T(2fj, + n + iX) : 2p: + n + iX = -m <^ RcA = 0, f3 = 2^i + n + m; ^ ' 

T(2fi + n + iv) : 2/i + n + iv = — m Rc^ = 0, a = 2/i + n + m . 

Here m > is an integer. From here on we always assume 0</3<a<2/i, so that only the first two 
lines in (|144[) remain possible. This assumption also prevents the occurence of the terms containing 
B n (v, A) in (|14ip . We now consider two special cases. 
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According to the above discussions, if N is a non-negative integer such that 2/i + 2N < a + (3 < 
2/i + 2N + 2, 

w ia (z, z')w ifj (z, z') = \ Kp iatil }{K)w K (z, z')dn 



R 



+ ^A n (ia, i{3)w l{2tl+2n _ a _ m . (145) 

71=0 

provided neither a nor (3 is an integer. If a + (3 < 2/i the formula holds without the A n terms. It is 
easy to check that Pia.ift^K*) _^ 0. For A n {i(y. 7 z/?) we find 

A n (ia, i(3) — 77i x 

n\2H^-R d - 2 Y (^-) r(-2/i - 2n + a + [3)T{-2n + a + f3)T(4fx + 2n-a-p) 

T(a - n)T(j3 - n)T(a + [3 - n)T(2fi ~a + n)r(2 j u - j3 + n)r(2 j u + n) x 
(-l) T T(n + 2/z-a-/3) 



T(2n + 2/i - a - (3) 



(146) 



All factors in this expression except the last fraction are positive since the arguments of the Gamma 
functions are positive because of the conditions 2/i > a > f3 and a + /3 — 2\i — 2n > 0. The last fraction 
is of the form 

(-l)T(7i + x) 



(-IT Hiq + X)- 1 . (147) 



T(2n + x) 

q 

The last product contains n negative factors and the result is positive. The positive sesquilinar 
form defined on SS(Xd) by the l.h.s. of (|145[) is thus, in this case, a positive superposition of 
positive sesquilinar forms, each corresponding to an irreducible unitary representation of the de Sitter 
group. 

It follows that a particle with parameter 17 (with < 7 < 2/i) can decay into two particles with 
parameters ia and i[3 (with the preceding conditions satisfied) provided 

7 = a + [3 - 2/i - 2n, i.e. 2/i - 7 = (2/i - a) + (2/i - (3) + 2n, (148) 

where n is a non-negative integer. 



C.2 Case v = ia, < a < 2/i, Im A = 0, Re A 7^ 

Only a particle in the principal series can decay into two particles with such parameters. There is no 
other restriction on the mass of the decaying particle. 
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